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We describe a first attempt to understand dynamical fermions within a Hamiltonian framework. As a testing 
ground we study compact QED3, which shares some important features of QCD4 such as confinement, glueballs, 
mesons, and chiral symmetry breaking. We discuss the methods used and show data for the chiral condensate. 



1. Introduction 

Nearly all lattice studies with dynamical 
fermions rely on pseudofermion techniques. This 
includes the standard Hybrid Monte Carlo 
method for non-local boson actions jjj and 
more recently multiboson methods proposed by 
Liischer O and extended to odd flavors by Borici 
and de Forcrand ||. While these methods are 
presently the most efficient available for dynam- 
ical fermion calculations, much improvement is 
still needed to combat algorithmic slowdown at 
light quark masses and large volumes. On the 
more theoretical side, one would also like to un- 
derstand the microscopic dynamics of quarks in 
lattice simulations. Ideally one hopes to see the 
effects of confinement on quark-antiquark pair 
separation and the effects of chiral symmetry 
breaking on light quark mobility. These ques- 
tions are most easily answered in a framework 
which does not integrate out the fermionic de- 
grees of freedom. With these goals in mind we 
investigate the possibility of simulating dynami- 
cal fermions in a confining gauge theory with ex- 
plicit fermion dynamics. 

In this proceedings article we describe the be- 
ginning of a research program at NC State to 
understand dynamical fermions within a Hamil- 
tonian framework. As a testing ground we 
study compact QED3, which exhibits confine- 
ment, glueballs, mesons, and chiral symmetry 
breaking. We discuss the methods that are be- 
ing used and show preliminary data for the chiral 
condensate and compare with strong and weak 
coupling results. 



2. Compact QED 3 

We follow the notation used in M and || . The 
Hamiltonian is 



where a is the lattice spacing and x the dimen- 
sionless strong coupling constant, 

x = e- 4 a- 2 . (2) 

W has three pieces, 

W = W + VxWi+xW 2 , (3) 

where 

l r 

(4) 

= £ »&(»*) [x t (r)U 1 (r) X (r + j) + h.c] 
f,j 

V 

and 771(7") = (-l) r2+1 , r} 2 (r) = 1. We use 
four component staggered fermions, which in the 
three-dimensional Hamiltonian formalism gives 
one flavor of fermion. The dimensionless mass 
parameter /1 is given by 

[A = 2me~ 2 . (5) 

3. Methods 

We start with a set of trial states |^-) and 
evaulate the quantities 

Mij = exp [-HT] , (6) 
M^=(* i |exp [-H(T+AT)}\^). 
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Improved trial states are given by 

|*;.)=exp[-!tfT]|^). (7) 
It is easy to see that 

M u = (*il^>> ( 8 ) 
M£ = (*;|exp[-#AT]|^}. (9) 

We now find a matrix T and unitary matrix U 
such that 

TijMj k T k i — da, (10) 
U i: jT^ k M^Ti m U mn = XiSi n . 

From these we can read off approximate eigenvec- 
tors 

\El)=Ti m U mj \*'j) (11) 
and eigenvalues 

\E l )=-^ T \ogX l . (12) 

The matrices M and M A are evaluated by di- 
viding into time slices, 

M ij = W | exp [-Mf] ■ ■ ■ exp [-Sf] |*$) . 

(13) 

Between the time slices we insert a complete set 
of states, and the basis we choose is the tensor 
product basis of gauge and fermion states. For 
the gauge states we use coherent states of the link 
fields, 

u j &® l> M(f)) = ^ {3) ® lr -\Mr)) 

(14) 

U}{x)® hr AMr)) = e-^^Mr)) . 

For the fermion states we use the usual position- 
space basis. For fixed coherent states of the 
gauge field between time slices, the fermion dy- 
namics are generated a classical background field, 

xHr)Uj(r)x(r + j) -» e iA ^*V(r)x(f + j). 

(15) 



4. Simulations 

In our calculations we keep the lowest few en- 
ergy fermionic states in computer memory and 
calculate their dynamics exactly by matrix mul- 
tiplication. We use a simple heat bath algorithm 
to propose new gauge updates. These updates 
are performed by sweeps through the spatial lat- 
tice, taking each time slice in succession. We use 
a Metropolis accept/reject decision based upon 
an estimator that includes the contribution of the 
low energy fermion states. In comparison with 
quenched simulations, the use of this estimator 
is slower by a constant factor proportional to the 
number of low energy fermion states. After each 
sweep through all time slices, the contribution of 
other fermionic states are calculated using diffu- 
sion Monte Carlo sampling. 

In our initial studies of small to mid-sized lat- 
tices we store in memory the nine lowest energy 
fermion states in memory, the strong coupling 
ground state and the eight translationally invari- 
ant states connected to the ground state by a sin- 
gle hop. In Figure 1 we show results for \i = 2. 
The data is extrapolated from spatial lattice sizes 
4 x 4, 8 x 8, and 10 x 10, and number of time 
slices 11, 15, and 19. The curves Sen s i> s 2, S3 
give the first few orders in the strong coupling 
expansion, and w is the leading behavior in the 
weak coupling expansion. We find behavior con- 
sistent with strong coupling estimates at small x 
and weak coupling estimate at large x. 

5. Conclusions 

We have described a first attempt to treat dy- 
namical fcrmions within a Hamiltonian frame- 
work. We considered compact QED3, which 
shares some important features of QCD4. At this 
time it is too early to say whether this approach 
would be viable for QCD4 with relatively light 
quarks. We are currently making accuracy and 
efficiency comparisons with pseudofermion meth- 
ods in QED3 H for exactly massless fermions. 

On the theoretical side there is clearly much po- 
tential in this approach to better understand the 
microscopic dynamics of quarks in lattice simula- 
tions. In future work we plan to study the effects 
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Figure 1. Chiral condensate for /i = 2. 
so.sij S2, S3 represent the lowest order estimates 
in the strong coupling expansion, and w is the 
leading result in the weak coupling expansion. 



of confinement and chiral symmetry breaking on 
quark dynamics. 

The author thanks Pieter Maris and Nathan 
Salwen for discussions. 
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